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Abstract 

Degenerate Hubbard models are studied using the Generalized-Gutzwiller- 
Approximation. It is found that the metal-insulator transition occurs at a 
finite correlation Uc when the average number of electrons per lattice site is 
an integer. The critical Uc depends sensitively on both the band degeneracy 
N and the filling x. A derivation is presented for the general expression of 
Ucix,N) which reproduces all previously known Gutzwiller solutions, includ- 
ing that of the Boson Hubbard model. Effects of the lattice structure on the 
metal-insulator transition and the effective mass are discussed. 

PACS numbers: 71.10.-hx,71.30.-hh,74.70.W 
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The simple Hubbard model, where each lattice site is occupied by at most two elec- 
trons, has been extensively studied in recent years as a model of strongly correlated electron 
systems.01i But in many systems, such as molecular solids Ceo and AajCeo,^ the conduc- 
tion band is degenerate which allows the occupation of more than two electrons per lattice 
site. For these systems, studying the degenerate Hubbard model on a lattice is the first 
step toward the understanding of electron-electronic correlations. Yet, there has been little 
progress since the pioneering works of Chao and Gutzwiller in early 1970s.i 

Recently, we succeeded in obtaining the Gutzwiller-like analytical solutions to degenerate 
Hubbard models using a scheme similar to the original Gutzwiller-Approximation.i A brief 
report of our results has been published.i The solution was used successfully to interpret the 
unusual metal-insulator transitions observed in fuUerides AajCeo-^S In this paper we present 
details of our derivation and the approximation. In addition, effects of the crystal lattice 
structure on the metal-insulator transition and the effective mass are also investigated. 

We consider the A^-fold degenerate Hubbard model defined on a lattice 

H= J2 t4,c^(^i,c. + ^ E ni^arii,^ , (1) 

<ij>=n.n.,a i,a^fl 

where < ij > are nearest neighbor lattice sites, a = (r, a) includes both the spin (a =1, |) 
and the orbital (r = 1, 2, ■ ■ ■ , A^) indices, and rii^a = cfaCi^a is the number operator. Both 
the Hubbard U and the nearest neighbor hopping t are assumed to be independent of lattice 
sites, electron spins, and molecular (atomic) orbitals. Let L be the total number of lattice 
sites, X the average number of particles per site. For simplicity, only total symmetric states 
are considered. Thus, the number of electrons per orbital per spin is m = xL/2N. For an 
integer x one expects an insulating ground state for f/ >> t. In such an insulating state 
there are x localized electrons at each lattice site. Hopping between sites are quenched due 
to the large correlation energy cost. On the other hand, for U = the system is metallic as 
electrons hop between lattice sites to lower the kinetic energy. One expects that there exists 
a critical Uc across which the system undergoes a metal-insulator transition. In general, the 
insulating ground state can be magnetically ordered. We will consider only non-magnetic 
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states. Though the quahtative results described below may be valid for magnetically ordered 
states, it is clear that further works are required to include the magnetic ordering. 

If an electron hops from a site with x particles to a neighboring site which also has x 
particles, then the correlation energy cost, U{Cl_i + Cj+i — 2C^) = f/, is independent of 
the band filling x and the degeneracy A^. The gain in kinetic energy is of the order of band 
width. Thus, in a simple mean field argument one might expect that Uc to be independent 
of the band filling x and the degeneracy A^. We will show that this is not true. In fact the 
solution we obtained in a Gutzwiller-type approximation predicts that Uc depends sensitively 
on both X and . 

Metal-insulator transition in non-degenerate Hubbard models have been extensively stud- 
ied in recent years.! At the half-filling (one electron per site) it is known: (1) In ID, the 
Bethe Ansatz exact solutionll shows that the ground state is insulating for [/ > 0. (2) In 
2D, the quantum Monte Carlo calculations suggest an insulating ground state for a mod- 
erate f/.i In addition, for an extended Hubbard model (the so called Kievelson-Schriffer- 
Su-Heeger model) it has been shown exactly that the Mott-Hubbard transition exists at a 
non-zero finite ?7c A and the Uc agrees with that obtained by the Gutzwiller- Approximation. 
(3) In 3D or higher it is generally believed that a finite Uc exists. (4) In the infinite di- 
mensions the Gutzwiller-Approximation is exact. From these results we conclude that the 
Gutzwiller-Approximation provides a reasonable estimate of the Mott-Hubbard transition 
in 3D degenerate Hubbard models. 

In the simple Hubbard model each lattice site are either empty, singly-occupied, or 
doubly-occupied. Within the Gutzwiller approach, one introduces a wave function which is 
a linear combination of states with different number of doubly-occupied sites with weighting 
factors to be determined variationally. To make analytical calculations possible, Gutzwiller 
made two approximations 1) In the thermodynamic limits, the expectation value of an 
operator in the variational function is approximated by the largest term contributing to it. 
This approximation is justified by the fact that the weighting factor, which includes various 
combinatorial factors, is a sharply peaked function. 2) In evaluating the expectation the 
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spatial correlations of doubly- occupied sites are neglected. This is equivalent to the mean field 
approximation. For the A^-fold degenerate Hubbard model, the problem is more complex as 
each site can be occupied by to 2N number of electrons. Without further approximation the 
combinatorial counting becomes intractable. However, close to the metal-insulator transition 
one expects that large fiuctuations from the average occupancy x is small as they cost large 
correlation energy. This leads us to make the additional approximation: 3) Each lattice 
site is restricted to three possible states: "empty" (E) — with x — 1 electrons, "singly- 
occupied" (S) — with X electrons, or "doubly- occupied" (D) — with x + 1 electrons. This 
approximation makes it possible to do variational calculations following similar steps of the 
original Gutzwiller- Approximation. We will refer to conditions 1), 2), 3) as the Generalized- 
Gutzwiller-Approximation (GGA) for the degenerate Hubbard model. 

Let r,fi,A be the total number of sites for types (D), (S), (E) respectively. Approxima- 
tion 3) and the conservation of particles require that L = r + fi + A, A = r. For a fixed 
r , there are many possible ways of participating L lattice sites. A given (D) site is distin- 
guished by identifying fiavors Pi,P2, ■ ■ ■ ,Px+i ^ (1? 2, 3, ... , 2N) of all x + 1 particles on the 
site. The numbers of lattice sites occupied by the same configuration in general depends on 
the configuration. For clarity, we will use the short-hand notation 7^ = 1, 2, 3, . . . , C2Pf^) 
for this number. There are €2^^^ = (a._)_i)!(2jv-z-i)! distinct possible configurations for each 
(D) sites. Similarly we define Ur and for sites of type (S) and (E). Respectively, there are 
Cjjv and C27v^ possible configurations for these sites. One has 

^ 7p = r, J2^r = n, E \ = A. (2) 

p=l r=l g=l 

Following Gutzwilleri, one introduces the weighting factor < t] < 1 for each (D) site, the 
GGA variational wavefunction {(pcGA > can be written as a linear combination of states 
with different {Xp}, {00,,}, {A^}, 

l^GGA >= E ivfMiK], Wr}, {AJ)|0{A,},V.},{A,} > • (3) 

{\p},{wr},{\q} 

The coefficient A is proportional to the total number of possible ways of participating the 
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lattice given {Xp}, {uJr}, {\}, 



^2N ^2N ^2N ^ 



\Ai{^p},{u;r},{X,})f ^ const X HCll x J] C^^r-E- 11 ^tr-n-^^ a. 

= const — — +1 —J , (4) 

where, the const is determined by the normahzation condition. Given the general wave- 
function Eq.(3), the expectation value of an operator is difficult to calculate. Within the 
Gutzwiller- Approximation one assumes that in the thermodynamical limit only the largest 
term contributes significantly to the expectation. For such an optimal term, the symmetry 
of the problem dictates that 



7p 


= 7 


p ir'x+i 








= A 


= NC^2N = 




(5) 




= CO 




(L-2r)/cf^. 





The relationship between 77 and F is determined by the largest term in < ■0gga|'0gga >■ 
This leads to 

^2(r+i) |^(^ + 1, 7, 7, • • • , 7; a; - 1, u; - 1, a;, • • • , a;; A + 1, A, A, • • • , A) |2 

= 77^^| A(7, 7, 7, • • • , 7; • • • , A, A, A, • • • , A) 1^ (6) 

Using Eq.(2), (4) and (5) one obtains 

^ v^A ^ F I {x + l){2N-x~+lj 
^ uj L-2T\ x{2N-x) ■ ^ ^ 

Keeping only the largest term in the normalization condition < ■0gga|'0gga >— 1 one finds 

(consty^ = r]^^\A{-f, 7, 7, • • • , 7; a;, a;, a;, • • • , a;; A, A, A, • • • , A)|^ (8) 

The optimal value of F (or rj) is determined variationally by minimizing the total energy. 
The energy consists of the kinetic and the correlation part. The kinetic energy is calcu- 
lated by evaluating the hopping term in Eq.(l). This requires the calculation of correlation 
functions 



]j,a) =< ^GGA\claCj^a\^GGA > ■ (9) 

Because each site can be occupied by x + 1, x, or x — 1 electrons only, the case i = j is easily 
evaluated to give 

p{t,a;t,a) = ^ 2N ' ^ ^ 

Where, in the last step Eq. (5) is used. For the case i ^ j there are three distinct possibilities: 
(a) None of two sites are doubly-occupied; (b) One of two is doubly-occupied; (c) Both are 
doubly-occupied. Summing all cases together leads to 



L{L - m) 

r(L-2r)x(2iv-x + i) / ^C2Vi y 

2NL^{2N-x) [ ^XC2N-ix-l) ^ ' 

r(L - 2r) 



x(2A^ - X + 1) + + l)(2iV - x 



2NL'^{2N - x) 

In deriving these results Eq.(5) and (7) are used. The kinetic energy per particle can be 
written as 

K.E. = Q{x,T,N)e{x), (12) 

where 

e(x) = ^ E <k) < (13) 

k<kF 

is the average kinetic energy per particle in the absence of the correlation. (The bare 
orbital energy has been chosen to be zero.) The quotient Q(x, r,A^),l which represents 
the reduction in kinetic energy of the correlated system in comparison with that of the 
non-correlated system, is given by 



Q(x,r,iV) 



T{L-2T) 



xL^{2N -X 

The correlation energy is determined by the expectation 



(^^Jx{2N -x + l) + ^{x + 1){2N -x)^ . (14) 



= L<^^T. (15) 

The same expectation when evaluated in the paramagnetic insulating state, where there are 
X electrons localized at each site, gives a value /, ^(^~^) , The difference, F, is the increase 
in the correlation energy in the metallic state. As the kinetic energy in the paramagnetic 
insulating state is zero, the increase in total energy per particle in the metallic state is, 

E{N, x) = Q{x, r, N)e{x) + , (16) 

xL 

Minimizing this energy with respect to F, one obtains the fraction of doubly-occupied sites, 
and the total energy per particle 

UJx,N)^ U 

Where, the critical U is given by 

^^x{2N-x + l) + ^{x + l){2N-x))\ 
Uc{N,x) = ^ ^|e»| . (19) 

As U increases toward Uc, the number of doubly-occupied sites approaches zero. For U > Uc 
the paramagnetic insulating state has a lower energy than the metallic state. Thus the GGA 
predicts a first order Mott-Hubbard metal-insulator transition at Uc- 

Eq.(19) is our main result. It shows that Uc depends sensitively on both the number of 
electrons per lattice site and the band degeneracy. In Fig.l we show several examples of 
this dependence with the band energy calculated from the nearest neighbor tight-binding 
Hamiltonian on the simple cubic lattice. A prominent feature is that Uc is the largest at the 
half filling, Uc{N, N) = 4(A^+l)e(A^). This implies that if the system is insulating at the half 
filling, then for all other fillings with integer number of electron per site the system is also 
insulating. This result has been used successfully to explain the unusual metallic/insulating 
properties of fuUeride materials A^^Ceo-^ 



The critical correlation also depends sensitively on the lattice structure through the band 
energy e{x). In Fig.2 we show examples of the phase diagram with e{x) calculated from the 
2-dimensional square, the simple-cubic, and the body-center-cubic lattice. As the number 
of nearest neighbor increases Uc decreases with respect to the band width. 

Our result, Eq.(19), reproduces all previously known analytical Gutzwiller solutions. 
The three cases that we are aware of are: (1) = 1, in this case the only commensurate 
filling is X = 1, Eq.(19) leads to the well known Brinkman-Rice criteria Uc = 8|e|.i (2) 
N = 2 and x = 1, this is the case studied by Chao and Gutzwilleri in early 1970s. The 
critical correlation they obtained, Uc = ^^^^y^|e|, is identical to that given by Eq.(19). 
(3) ^ oo with a finite x, in this case the fermion problem is equivalent to the boson 
problem. The boson Hubbard model has been widely studied in recent years0^@. Our 
result, Uc = {y/x + v^ar+T)^|e(x)|, is the same as that obtained for the boson Hubbard 
model using the Gutzwiller- Approximation.0 It is also close to the value given by quantum 
Monte Carlo calculations.El This suggest that in the limit of high band degeneracy our 
Generalized-Gutzwiller-Approximation is quite accurate. It would be very interesting to do 
Monte Carlo calculations explicitly for the degenerate model and compare with our results. 

The physical properties of the system can undergoes dramatic changes when the metal- 
insulator transition is approached. These include the magnetic susceptibility and the trans- 
port effective mass. As Brinkman-Rice first pointed out,i the effective mass m* is critically 
enhanced near the metal- insulator transition. The enhancement is proportional to 1/Q. 
Using Eq.(14) and (17) one finds 

rrf^]_^ S\e{x)\ 

m, Q xUc{x,N){l-{U/UcY)' ^ > 

where is the effective band mass. One observes that the effective mass diverges quadrat- 
ically as the metal-insulator transition is approached. However, the enhancement is less 
dramatic for large A^ or large x. Brinkman-Rice also showed that the magnetic susceptibil- 
ity also diverges as 1/(1 — (f/ /f/c)^). We suspect that similar result holds for the degenerate 
Hubbard model. However, this remains to be proven as our calculation is valid only for 
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non-magnetic states. 

In conclusion we present results of analytical studies on degenerate Hubbard models 
using the Generalized- Gutz wilier- Approximation. It is shown that for any filling with integer 
electron per lattice site there exists a critical correlation energy above which the system is a 
Mott-Hubbard insulator. The general expression we found for Uc depends sensitively on the 
band degeneracy, the number of electron per site and the lattice structure. It also reproduces 
all previous know Gutzwiller solutions as special limits. 
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FIGURES 
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FIG. 1. Phase diagrams of metal-insulator transitions in degenerate Hubbard models as 



predicted by Eq.(19). Critical correlation Uc (scaled by the band width) is plotted against the 



band filling x (scaled by 2N). The band energy is calculated from the n.n. tight-binding model 



on the simple-cubic lattice. Shown are results for A'' = 1,2,3,4,5. The Uc is the largest at the 



half-filling x = N. In both Fig.l and Fig.2 only points, corresponding to integer number of electrons 



per site, are meaningful. Lines are draw using Eq.(19). 
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FIG. 2. Phase diagrams of metal-insulator transitions in degenerate Hubbard models for 
different lattice structures. Shown are three examples with N = 5 and the band structure calculated 
for the 2-D square, 3-D simple-cubic, and 3-D body-center-cubic lattices with the nearest neighbor 
tight-binding hamiltonian. Note that given the same band filling and the degeneracy, the critical 
correlation can changes substantially from one type of lattice structure to another. 
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